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Deflection angles of massive test particles moving along an unbound trajectory in the 
Schwarzschild metric are considered for the case of large deflection. We analytically consider the 
strong deflection limit, which is opposite to the commonly applied small deflection approximation 
and corresponds to the situation when a massive particle moves from infinity, makes several revolu¬ 
tions around a central object and goes to infinity. For this purpose we rewrite an integral expression 
for the deflection angle as an explicit function of the parameters determining the trajectory and 
expand it. Remarkably, in the limiting case of strong deflection, we succeed in deriving for the first 
time the analytical formulas for deflection angles as explicit functions of parameters at infinity. In 
particular, we show that in this case the deflection angle can be calculated as an explicit function 
of the impact parameter and velocity at infinity beyond the usual assumption of small deflection. 
PACS numbers: 04.20.-q - 04.20.Cv - 04.25.Nx - 04.20.Jb - 95.10.Ce - 95.10.Eg 


I. INTRODUCTION 

The motion of massive test particles in the 
Schwarzschild metric has been extensively studied 
throughout the years [Ml- Depending on initial param¬ 
eters, there are different types of trajectories possible. 

Let us consider an unbound orbit with a massive par¬ 
ticle moving from infinity to a central object and then 
again to infinity. If the impact parameter of the particle 
is large, the trajectory is an almost straight line with a 
deflection by the small angle a <C 1 (see, for example, 
Ref. @): 

d = ^(l + 2^, R S = 2M, G = c= 1 . ( 1 ) 

Here Rs is the Schwarzschild radius, and b and v are the 
impact parameter and velocity of incident test particles, 
respectively. 

Beyond this simple case, the deflection angle of a mas¬ 
sive particle moving from infinity to a central object and 
then to infinity cannot be calculated in such an easy way. 
The exact formula for the deflection angle is written as 
an integral from R (the distance of the closest approach) 
to infinity and can be expressed via elliptic integrals. 

In this paper we pay our attention to another limit¬ 
ing case in which the deflection angle can be calculated 
analytically, the strong deflection limit (d 1), which 
corresponds to the situation in which a massive particle 
moves from infinity, makes several revolutions around a 
central object and goes to infinity. The deflection angle 
goes to infinity when the distance of the closest approach 
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tends to the radius of unstable circular orbit of the parti¬ 
cle (which depends on initial parameters). Thus, a strong 
deflection limit implies that the parameters of the orbit 
are close to their values at this radius. 

Analytical formulas for strong deflection limit were de¬ 
rived only for the case of photon motion [§]. The deflec¬ 
tion angle diverges logarithmically when R approaches 
3 M. For photon deflection beyond the small deflection 
limit, see classical papers and books i, 1 , i, 0 , 1 , 
a a- The subject of photon deflection in a strong de¬ 
flection situation is very popular in gravitational lensing 
theory. Photons that perform several revolutions around 
the central object can form so-called relativistic images 
of source. For gravitational len sing in the strong deflec¬ 


tion case, see Refs. [D 
EOl, (H, m, and m 


m, 0 , 0 , 0 , 


The integral in the exact formula for the deflection an¬ 
gle of massive particles includes a polynomial of the third 
order in the radicand, and the expression via elliptic in¬ 
tegrals is usually being written using the roots of this 
polynomial. In the general case, these three roots are 
different and supposed to be found numerically for given 
external parameters E (energy at infinity per unit rest 
mass) and L (angular momentum per unit rest mass), 
which determine the trajectory and R. For given E and 
L, one has to obtain the roots of the polynomial numer¬ 
ically before being able to calculate the deflection angle 
(see, for example, Refs. [Ij and Q). Thus, if the im¬ 
pact parameter is not large, the deflection angle is not 
expressed explicitly via external parameters E and L. 

The divergence of the deflection angle leads to diffi¬ 
culties in the derivation of asymptotic formulas for large 
angles, and it seems that the most appropriate way is 
to expand the exact integral expression. For this pur¬ 
pose we need the integral to be an explicit function of 
parameters determining the trajectory. To achieve it, we 
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suggest using either the pair (R, L) or (R, E ) as an in¬ 
dependent pair of parameters, instead of (E, L). In this 
case, all roots can be found analytically, and therefore 
we can perform all the calculations analytically. In this 
paper we obtain the deflection angle in the form of ellip¬ 
tic integrals as an explicit function of the two parameters 
( R , L) or ( R , E ). A similar way is used in the consider¬ 
ation of photon deflection, in which a similar polynomial 
is expressed via R. 

Expanding the integrals in the asymptotic case of large 
deflection angles, we obtain the formulas for the deflec¬ 
tion angles as analytic functions of (R, L) or ( R , E). 
Nevertheless, for applications it is more convenient to 
determine the trajectory by the parameters at infinity. 
The distance of the closest approach R is not a parame¬ 
ter at infinity, and therefore it is not good as an external 
parameter. Luckily, in the strong deflection limit, it be¬ 
comes possible to express the deflection angle in terms 
of the impact parameter b instead of R. Thus, only the 
parameters at infinity are used for calculation of the de¬ 
flection angle, which is convenient for applications. A 
similar approach was used in the description of photon 
deflection in Ref. [HI . 

In this work we present for the first time the analytical 
formulas for the deflection angles of massive particles in 
the strong deflection limit (a 1) as an explicit function 
of external parameters at infinity, (6, L) or ( b , E). Pa¬ 
rameters at infinity are simply related with each other, 
so it is possible to calculate the deflection angle with our 
formulas in very different situations. For example, our 
formulas allow one to calculate analytically the deflec¬ 
tion angle of a particle with a given impact parameter 
and velocity at infinity as initial conditions, beyond the 
usual small deflection case. 

Results of the present paper can be useful for the inves¬ 
tigation of star motion around supermassive black holes. 
In this case, the ratio of masses is extremely high, and 
such orbits now are often referred to as extreme mass ra¬ 
tio inspirals (EMRIs). Interest in these orbits is very high 
now, since they are some of the most prominent sources 
of gravitational waves (see, for example, the pioneering 
work of Zel’dovich and Novikov H). With appropriate 
initial conditions, the test body can perform several rev¬ 
olutions around the central object. In this case, weak 
field approximation is not valid, and the usual way is to 
solve the strong field equations numerically. Our formu¬ 
las allow an easy explicit calculation of the angular de¬ 
flection of a star or compact object in EMRIs and then 
finding, for example, the number of revolutions of the 
test body around a supermassive black hole for various 
initial conditions, which can be useful for gravitational 
wave estimation [ 25j . 

In our recent work |24j , we considered unbound photon 
orbits in the Schwarzschild metric in presence of homoge¬ 
nously distributed plasma, in the case of strong deflec¬ 
tion. In particular, we have derived the formulas for the 
deflection angles, in the case of strong deflection. In the 
paper of Kulsrud and Loeb [26| , it was shown that in ho¬ 


mogeneous plasma the photon wave packet moves like a 
massive particle with velocity equal to the group velocity 
of the wave packet, the rest mass equal to the plasma fre¬ 
quency, and its energy equal to the photon energy. This 
analogy allows one to apply the results of that paper 
to the calculation of deflection angle of massive particle, 
with a corresponding change of variables. In some places 
of the present paper (where possible), we will refer to 
some results of that paper. For the subject ofgravita- 
tional lensing in plasma, see also Refs. [27], [28], [29j . 
01, 0, 01, and Up. 

The paper is organized as follows. In Sec. II we derive 
the well-known exact expression for the deflection angle 
of the massive test particle moving in the Schwarzschild 
metric along an unbound orbit. In Sec.Ill we express 
the formula for deflection via elliptic integrals, working 
first with the pair of independent parameters (i?, L ) and 
then, separately, with the pair (R, E). In Sec. IV we find 
a critical value for the closest approach distance at which 
a massive particle going from infinity remains on the cir¬ 
cular orbit with an infinite number of circles around the 
center. This critical value depends on L (or on E). The 
deflection angle goes to infinity when the distance of the 
closest approach goes to this value. In Sec. V we derive 
an analytical formula for the particle deflection angle in 
the Schwarzschild metric, in the limit of strong deflection. 
In Sec. VI we conclude and discuss our results. 


II. EXACT DEFLECTION ANGLE 


Let us consider a massive test particle moving in the 
Schwarzschild metric: 

ds 2 = —A{r ) dt 2 + f \ + r 2 (dO 2 + sin 2 9dip 2 ) , (2) 

A(r) 


A(r) = 1-. 

r 

The system of units used here is 
G = c = 1, the Schwarzschild radius Rs = 2 M. (3) 
Indices are 

i,k = 0,1,2,3; a,/3 = 1,2,3 (r,0,tp), (4) 


signature {-,+,+,+}. 


( 5 ) 


Equations governing the orbit of the test massive par¬ 
ticle (see, for example, Ref. Q) are: 

( — J =E 2 — V 2 (r) (radial part of motion), (6) 


dip L 
dr r 2 


(angular part of motion). 


( 7 ) 
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Here E is the energy at infinity per unit rest mass of the 
particle, L is the angular momentum per unit rest mass 
of the particle, and the effective potential V (r) is 


H 2 (r) = 




( 8 ) 


III. EXPRESSION OF THE DEFLECTION 
ANGLE VIA ELLIPTIC INTEGRALS 

Here we express the integral for the angle m via el¬ 
liptic integrals. 

Let us introduce notations 


From equations of orbit, we can obtain the equation of 
the trajectory as 


dip L 1 

*“ ? V /E s -A(r )(! + #)' 


(9) 


— = u, — = uq , A(u ) = 1 — 2Mu, A(uq) = 1 — 2Muo, 
r K 

(13) 

and rewrite Eq. m in terms of u instead of r. We will 
refer to the expression in radicand in Eq. m as /(it). 
Equation CD becomes 


Assume that the test particle moves in such a way that 
its p coordinate increases. Then the plus sign in Eq. 
© corresponds to motion with the coordinate r also in¬ 
creasing, and the minus sign corresponds to motion with 
decreasing r. 

For a particle which moves from infinity to the dis¬ 
tance of the closest approach R (minimal value of the 
coordinate r) and then to infinity, a change of angular 
coordinate is 


Aip 



dr 

\]e 2 - A(r) (1 + £) 


( 10 ) 


u 0 


a = 2 


Ldu 

\/7M 


(14) 


where f(u) = E 2 — (1 — 2Mu)(l + L 2 u 2 ). (15) 

Coefficients of the polynomial f(u) are functions of E and 
L , and the roots of this polynomial could be found either 
numerically or by a very complicated analytical solution 
of the third-order algebraic equation. We suggest using 
R and L (or R and E) as another pair of independent 
parameters. Using the boundary condition m , we write 

E 2 = (1 — 2Mu 0 )(l + L 2 Uq) , (16) 


+ 



dr 

yj E2 - A(r) (1 + £) 


= 2 


OO 


R 


L 

^2 


dr 

yj E 2 - A(r) (1 + £) 


The motion along the straight line corresponds to the 
change of the angular coordinate A <p = tv. Therefore the 
deflection angle can be written as 


a = 2 


dr 


\!E 2 - A(r) (1 + £) 


— 7r. 


(id 


The closest approach distance R and the parameters E 
and L are connected at the point r = R by the following 
boundary condition: 


/(it) = (l-2Mtto)(l+T 2 Uo)-(l-2Mit)(l+TV). (17) 

Now the coefficients of the polynomial are functions of R 
and L , and one of the roots u = 1/R is evidently visible. 
The other two roots (as functions of R and L) can be 
found easily by solving the quadratic equation. This ap¬ 
proach is productive because we can perform all further 
calculations analytically and obtain the deflection angle 
in the form of elliptic integrals as an explicit function of 
the pair of parameters R and L. A similar approach is 
usually applied for the deflection of photons in vacuum, 
where an analogous polynomial is expressed via R (see, 
for example, Ref. [17j | and other references in the present 
paper). We would like to mention that in the case that 
is of relevance for the subject of the paper the polyno¬ 
mial /(it) has three real roots, and R is the biggest of the 
three values for which f(\/R) = 0. Finally, 

/(it) = 2 ML 2 [u - ua){u - u B ){u - uc), (18) 


E 2 = A(R) (l + 1 (12) 

The trajectory and the deflection angle are thus com¬ 
pletely determined by any two parameters from { R , E, 
L }, with the third parameter being expressed through 
Eq. (HU). To compare, the photon motion is determined 
by only one parameter, either R or b , which are uniquely 
connected with each other. 


where 

R-2M + Q 1 R-2M-Q 

UA ~ 4 MR ’ UB ~ U ° ~ R ’ UC ~ AMR 

(19) 

e 2 = - 2M > 2 + SMIR - 2M > (‘ - L ™-m ) ■ 

( 20 ) 
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For the deflection angle d, we obtain the expression 

u 0 

2 f du 

V2 M J y/{u — ua){u - ub)(u - uc) 

Using Ref. [34J, we write 


( 21 ) 


and y given by Eq. (1251) satisfy the relation (EZD, inde¬ 
pendently on the form of Q. Therefore, we obtain the 
deflection angle d in the form 


d = 4 



[F(l,fc) - F(z, k)\ - 7r. 


(29) 


du 


y/(u — Ua){u — Ub){u — Uc) V U A — uc 


F{z,k ) 


( 22 ) 


_ / (UA ~ Uq)ub , _ j uB - UC 

V ( U B ~ U C )u A ’ V U A -U C 


[UA > UB > 0 > Uc] • 

Here F(z, k) is an elliptic integral of the first kind [34], 

v 


F( ‘- k) = Ijn 


dx 


dO 


o \/(l - x 2 ){l - k 2 x 2 ) J y/i -k 2 sin 2 0 

(23) 

= F(ip, k) , F(sin ip, k) = F(ip,k ), x = sin 0 , z = sin <p . 

Using Eq. E21) in Eq. (12T1) and substituting Eq. (fT!?l) . 
we obtain the deflection angle a in the form 


/ R 

ot = ^\l-Q F (y,k)-i r, 


(24) 


where y = 


8 MQ 


(6 M -R+ Q)(R - 2 M + Q) 


. (25) 


k = 


IQM — R + Q 


2 Q 


(26) 


The expression (l24l) can be also written differently, due 
to the property of elliptic integrals [36], which in our 
notations reads 


F(z,k) + F(y,k) = F(l,k), (27) 


provided y/l — k 2 = ^ = 1. 

Vl - \/l — y 2 

Here E(l, fc) = F(n/2, k ) = A'(fc) is the complete elliptic 
integral of the first kind. It is easy to check that 


Thus, we have obtained, for the first time, the formulas 
for the deflection angle m , m of massive test parti¬ 
cles, where arguments of elliptical integrals are expressed 
explicitly via parameters R and L, defining the trajec¬ 
tory. 

To obtain the deflection angle as a function of R and 
E , we rewrite Eq. m as 

L2 =* 2 {im - 0 • (30) 

and substitute Eq. (1301) into the expression for the de¬ 
flection angle m- The deflection angle then will be in 
the form of Eq. Cl with the expression f(u) in the form 


^ 2Mi? 2 (RE 2 — R + 2M) 

JW = - R-9M - {u-ua){u-ub){u-uc) ■ 


where 


R-2M + Q 1 

UA= AMR ’ UB = U0= fl lUC = 


( 31 ) 

R-2M-Q 


AMR 


(32) 


Q 2 = {R- 2 M) 2 + 8 M(R - 2 M) 


1 


1 , 2 M -• ( 33 ) 

1 + R(E*~ 1 ) 


Finally, in this case we obtain the same expressions 
and (EH1) with the only difference in the definition of Q. 
For more details about this case, see our work in Ref. 


We obtain the exact expression for the deflection angle 
as an explicit function of the pair ( R,L ) or ( R,E ), see 
Eqs. (1241) and m, with the variables y given by Eq. 
K51) . k by Eq. ESI) , z by Eq. ESI) , and Q in the form of 
Eq. EU1) for the (R,L) case or Eq. E3l) for the ( R,E ) 
case. 

The exact expression allows one to calculate the deflec¬ 
tion angle for any ranges of a. In case of small deflection 
angles (d <C 1), there is an analytical formula for the de¬ 
flection angle CD- In the opposite case of strong deflec¬ 
tion (a> 1), we can also obtain an asymptotic analytical 
formula for the deflection angle, expanding formula ESI) 
for the exact deflection. 


IV. CRITICAL DISTANCE OF THE CLOSEST 
APPROACH FOR UNBOUND MASSIVE TEST 
PARTICLES 


2M + Q-R 
6 M+Q-R 


( 28 ) 


To obtain an expression for the angle in the limit of 
strong deflection, we need to expand the integral ESI) 
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about the value of R , at which the deflection angle goes 
to infinity. It will be referred to as the critical value of 
R. 

The deflection angle q can be written in terms of an 
effective potential 0, 


a = 2 


Ldu 


/E 2 - V 2 (u) 


where 


V\u) = (l-2Mu)(l + LV). 


From the condition 


dV 2 (u) 
du 


= 0 


we obtain 


u = 


1 ± s/1 - 12 M 2 /L 2 


6 M 


(34) 

(35) 

(36) 

(37) 


The maximum and minimum of the effective potential 
are located at the points Tm and r m , respectively: 


tm = 


6 M 
1 + x 


6 M 
1 — x 


. 12 M 2 

x= \l 1- y /~) r M < r„ 


(38) 

The point r = tm corresponds to an unstable circular 
orbit of the particle, and r m corresponds to a stable cir¬ 
cular orbit. The critical distance of the closest approach 
(at given L) is being reached when the value of R coin¬ 
cides with tm- In other words, the unbound orbit with 
infinite deflection takes place if R = tm- For simplicity, 
we will denote the critical value of the closest approach 
as tm- 

The critical distance of the closest approach tm is a 
function of L. To derive the strong deflection limit as 
a function of (R,L), one should expand the integrals in 
exact formulas for deflection angles near R = vm where 
tm is expressed in terms of L. Analogously, to get the 
strong deflection limit as a function of (R,E ),one needs 
to rewrite tm as a function of E instead of L. 

Using Eqs. (l30l) and (1551) . we express the critical dis¬ 
tance tm for unbound massive particles coming from in¬ 
finity as a function of E in the form 


3E 2 — 4 + 3 E 2 y 
™ = 2(g-l) (39) 

Noticing that 3 E 2 — 4 + 3 E 2 y = ^E 2 (3y — 1)(1 + y), we 
write tm as 


3 (3y — 1)(1 + y) ,, 
rM = 4 —fm— M ' 

1 E 2 


(40) 


Then, because of the equality E 2 (3y—l)(3y+l) = 8 (U 2 — 
1), tm takes the following form: 


r M = 6M 


1 +y 
1 + 3 y 


y = 


= a/1 — 


8 

9 E 2 


(41) 


At given L (or E ), the trajectory is determined by the 
choice of the second parameter R. The deflection angle 
goes to infinity when R approaches tm, and the massive 
particle performs an infinite number of turns at the ra¬ 
dius tm- Therefore, is the critical (minimal) distance 
of the closest approach R , for given L (or E). For un¬ 
bound orbits of massive particles coming from infinity, 
making several loops around a black hole and going back 
to infinity, the value of the critical distance of the closest 
approach is between 3 M and 4 M. 

We should mention here that the points of maximum 
tm and minimum r m of the effective potential appear 
only for L > 2\[3M. According to Eqs. (1551) and (I4T1) . 
tm varies from 6 M to 3 M if L varies from 2s/3M to 
infinity or, equivalently, if E varies from -s/%/3 to infin¬ 
ity. In this paper we consider only the case in which the 
particle moves from infinity; therefore, for the studied 
orbit, it has E > 1. As we are interested in the case in 
which the particle reaches the point R and then returns 
to infinity, we need V/m) to be more than unity. In 
another words, the point of the maximum of the effec¬ 
tive potential should be higher than the value of effective 
potential at infinity. This happens if L > 4M (see Ref. 
[8|). Therefore, for the studied orbit, r« varies from 4M 
to 3 M if L varies from 4 M to infinity or, equivalently, if 
E varies from 1 to infinity. 


V. DEFLECTION ANGLE IN THE STRONG 
DEFLECTION LIMIT 

A. As a function of R and L 

When R tends to tm ( R > Ovt), the deflection angle 
goes to infinity. To obtain the analytic expression for the 
deflection angle in the strong deflection limit (d 1), 
one needs to expand the integral (l29l) around the point 
R = r M - 

As follows from Eqs. (l26l) and (l20l) . k ~ 1 corresponds 
to R ~ 7'm , so we can use the expansions of the elliptic 
integrals about k = 1 [HI, +/]: 


F{ 1, k) — In 


VT^k ^ ’ 


F(z, k ) ~ In tan 


arcsin(z) tt 


(42) 


(43) 


The expression 
to the following: 


2 4 / 

for the deflection angle is simplified 


(44) 


[r 

In 4 

m 1 + * 1 

/ Q 

VI -k 2 



Expanding (1 — k 2 ) near R = tm /m is defined by Eq. 
@], after simplifications (see Appendix A) we have 


9 Mx 


(45) 
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Thus, expanding Eq. (l44l) near R = tm and keeping the 
leading terms, we obtain 


a = - -= In 

\Jx 


(R - r M )M (1 + z) 


21 


4r 2 M x 


1-z 2 


(46) 


where z [expression (1281) with Q taken in the form of Eq. 
m] is simplified to 


2 = 


2x- 1 


3x 


(47) 


Substituting z, we obtain the deflection angle a as a func¬ 
tion of the closest approach distance R and L in the form 


a(R , L) =- -= In 

Va; 


Zi (x) 


R~r M 


tm 


(48) 


5a; - 1 + 2\j3x(2x - 1) 

a <*> = - ui -■ (49) 


where 


and tm and x are defined in Eq. (1881) . This formula 
is asymptotic and valid for R close to r«- 


B. As a function of R and E 

In this case, it follows from Eqs. (l26l) and (fTTTTl) that k ~ 
1 corresponds to R ~ tm , so we use the same formulas 
m, 63 > 63- Expanding (1 — k 2 ) near R = xm [xm is 
give n by 6D]> after significant simplifications (see Ref. 
[24 ] for details) we have 


1 - k 2 ~ 2 


1 +y M(R-r M ) 


y 


(50) 


' M 


Expanding Eq. (PHI) near the point R = xm and keeping 
the leading terms, we obtain 


a = — 2\ 


/ 1 + y i 

"1 + 2/ ( R-x m )M (1 + z) 2 ' 

V 2 , 

l 8y x 2 m 1 — 2 2 J 


- 7T, 


where z [expression 
is simplified to 


(5i) 

with Q taken in the form of Eq. 


C. As a function of b and L 

The impact parameter b for a massive particle is de¬ 
fined uniquely by E and L as f|J, [g] 


b 1 = 


E 2 - 1 ' 


(55) 


Let us assume that b cr is the minimal impact param¬ 
eter, corresponding to the critical distance of the closest 
approach 7'm . Since we need b to be a function of R 
and L, we substitute E 2 from Eq. (fl^l) into Eq. (l55l) . 
Expanding b in Eq. (1551) around the point R = tm, w e 
obtain 


b = b cr (x) + bi(R — tm) 2 


where 


b cr (x) = r M 


2x- 1 


h = 


3x 


2(2x 1 )r M 


b C r ( X ) . 


Writing 


(R - x M ) 2 _ 2_2 b-b cr (x) 
r 2 M Z b cr {x) 


(56) 

(57) 

(58) 

(59) 


we obtain the deflection angle a as a function of the im¬ 
pact parameter b and L, 


a(b, L ) =- -f= In 

\Jx 


2 z 2 z 2 (x) 


b- b cr (x) 
bcr (jx) 


(60) 


where z is given by Eq. (S3, Zi(x) is given by Eq. 63), 
b cr {x) is given by Eq. (1571) . and tm and x are defined 
in Eq. (1551) . This formula is valid for b close to b cr (x), 
where b cr (x) is the critical value of impact parameter at 
given L (1571) . 


D. As a function of b and E 


z = 


3y-l 

6 y 


(52) 


Substituting z, we obtain the deflection angle a as a func¬ 
tion of the closest approach distance R and E in the form 


a(R,E) = -2 
where 


1 + y 
2 y 


In 


*2 (y) 


R-x m 


XM 


(53) 


9 »- 1 + VM j izA, (54) 

48 y 


and rM and y are defined in Eq. (PHI) . This formula (l53l) 
is asymptotic and valid for R close to tm- 


In this case we need b to be a function of R and E, and 
we substitute L 2 from Eq. (15U1) into Eq. (1551) . Expanding 
b in Eq. (1551) around the point R = r«, we obtain (for 
this case, see Ref. [24]]) 

b = b cr (y) + h(R- r M ) 2 , (61) 

where 

b C r(y) = V3 x m \J i ( 62 ) 


3-\/3 y I 1 + y 
2 x M V 3y - 1 


3 y 

2 r M 


bcriy) ■ 


(63) 
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Writing 

(R - r M ) 2 _ _2_ b — b cr (y) 

r M 3 y b cr (y) ’ 1 j 

we obtain the deflection angle d as a function of the im¬ 
pact parameter b and E: 


a (b, E) 


J 1 + y in 

"2 z 2 {y) b-b cr {y)' 

V 2y 

3 y b cr (y) 


(65) 


where z%(y) is given by Eq. (151[1 . b cr (y) is given by Eq. 
(ED, and tm and y are defined in Eq. (ED- This formula 
is valid for b close to b cr (y ), where b cr (y) is the critical 
value of impact parameter at given E (1621) . 

If impact parameter b and asymptotic velocity at infin¬ 
ity v are used as initial parameters, the deflection angle 
can be calculated by the formula (1651) with 


E = 


1 


\/l — v 2 


In other cases connection L = Evb (see, for example, Ref. 
H) can be used. 


formula one should calculate the new pair of variables 
using the exact relations, for example, Eq. ED- 

In the limiting case of L —> oo or E —> oo, our formulas 
for massive particles transform to formulas for photons in 
vacuum; for the deflection angle of photons in the strong 
deflection limit with R, see Ref. [J|, and for the formula 
with b , see Ref. [l5| . 


VI. DISCUSSION 

The main result of present paper is a derivation of the 
analytical formulas (1551) . ED, ED, lt65l) for deflection 
angles of a massive test particle in the strong deflection 
limit (a> 1). The large deflection limit is the opposite 
of the commonly used small deflection case and works 
very well for the situation in which a massive particle 
goes from infinity, performs several revolutions around 
the central object ,and then goes to infinity. Formulas are 
written as functions of parameters at infinity determining 
the trajectory. 


E. Comparison of formulas 

We have shown the derivation of the formulas for the 
two cases {(R, L), (6, L )} and {(R, E ), (6, E)} from the 
very outset to the final results. In the first case, we used 
L as the second independent variable, and in the other 
case we used E instead. 

Here we would like to discuss also how these two cases 
can be transformed to each other by the appropriate 
change of variables. 

In the case of formulas with the distance of the clos¬ 
est approach R, Eqs. ED and ED, ^ is easy to do. 
Comparing Eqs. (1551) and ED- one easily finds that 
x = 2j//(l + y). Substituting this into Eq. (1551) . we 
obtain the deflection angle d as a function of the closest 
approach distance R and E in the form (1551) . 

In the case of formulas with the impact parameter b , 
Eqs. (15U1) and ED, the change of variables is more com¬ 
plicated and cannot be performed by this simple substi¬ 
tution. In the derivation of these formulas, we expand 
the impact parameter up to quadratic terms (R — tm) 2 - 
Therefore, in this case it is necessary to perform the 
change of variables with the same accuracy. To get the 
expression ED from Eq. ED, it is necessary to sub¬ 
stitute the exact expression of L via R and E [see Eq. 
ED] to Eq. (1561) and expand it up to quadratic terms 
{R — xm) 2 - In this procedure, additional quadratic terms 
will arise from the term b cr {x ) in Eq. ED, and after 
simplifications we can finally get the expression ED- it 
is also possible to get the expression (1551) from Eq. ED- 
Therefore, we can transform formulas (1601) and (1651) into 
each other. 

We should emphasize therefore that if one wants to 
use other variables instead of the pair used in the given 
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Appendix A: Calculations for Sec. V 

To find an expansion ED of (1 ~ fc 2 ) near the point 
R = tm, let us rewrite k 2 as 

, 2 6 M-R + Q 1 16 M-R , A . 

k = -2Q-=2 + 2 —Q— (A1) 

For convenience, let us calculate Q 2 /(6M — R) 2 . Here 
Q 2 is 

« 2 = <*- 2 M > 2 + -m(i- L a™ R l M) ) - 

(A2) 

Let us write R as 


R — f’M H - AR . 


(A3) 















Here A R = R — tm is a small variable: 


(A4) 


Expanding (6 M — R) 2 near R = vm (small variable here 
is A R), we obtain: 


A R <C tm , 


and 


tm = 


6 M 
1 + x 


(AS) 


Expanding Q 2 near R = r« (small variable here is 
A R), we obtain: 


where 


Q 2 — Q o + QiAR , 


3 6M 2 x 2 4Mx(1 + 4x) 

Q 0 = 77 —-77, Qi = 


(6Af - RY ~ 


2 36M 2 x 2 (1 + x 


(1 + a;) 2 \ 3Mx 


1-^TTTAi? . (A6) 


Keeping the leadings terms, we obtain Q 2 /(6M — R) 2 in 
the form 

Q 2 4(1 + x) 2 

t---^~1 + —-—A R. (A7) 

(6 M - R) 2 9 Mx K ’ 


With using (ED we obtain finally 


1 - k 2 ~ 


(i + xY 

9Mx 


A R. 


(A8) 


(l+z ) 2 


1 + x 
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